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In the past year, we have made considerable progress in a number of areas including algorithm devel- 
opment. completion of two major case studies, and the development of a new EUV flux model. As a result, 
there have been a major improvement in our ability to model global emissions in support of NASA’s imaging 
plans. Five papers are current y in press. One of these was written by an undergraduate student. 


£r' 


1. HIGHLIGHTS 

• Developed a new algorithm to allow physical models to reproduce observed NmF 2. 

• Investigated the relationship between NmF2 and FI0.7 at Millstone Hill during 1990. 

• Developed a new solar EUV flux model. 

• Statistical survey of anomalously high nighttime electron T e at Millstone Hill. 

• Conducted a case study of the March 1990 magnetic storm. 

• Conducted a comparison between theory and data of magnetically quiet behavior of the winter iono- 
sphere at Millstone Hill. 


2. DESCRIPTION 

The ability to accurately model F-region electron densities is extremely important in all aspects of 
ionospheric modeling. Our FLIP model has proved capable of reproducing the measured daytime densities 
very well using a technique that we developed to eliminate uncertainties due to neutral winds [ Richards , 

1991]. However, in common with other models, the nighttime densities are generally about a factor of 3 
smaller than the observed densities. 

There are many applications for which the electron density is just an ancillary input to the calculation 
to some other quantity to be studied, for example, when calculating emission rates or electron temperatures. 

These applications require ve - y accurate electron densities. Empirical models such as the IRI model repro- 
duce the average measured N mF2 and hmF2 very well but are unable to produce day-to-day variability the 
way a physical model can. Also, the topside density profiles are poorly modeled mainly because the topside 
is poorly sampled and more variable. 

We have developed a new technique which allows the FLIP model to reproduce the measured NmF 2 
as well as the measured hmF'l (or an empirical model NmF2 and hmF2). The FLIP model then fills in the 
topside ionosphere self-consistently. Thus the accuracy of all those densities, temperatures, and emissions 
that depend on electron density has been greatly improved. The contour plot in Figure 1 demonstrates the 
new capability of the FLIP model to reproduce the observed electron density. A paper on this new technique 
is in preparation. The IRI model has been modified to make it a callable routine that returns just NmF 2 and 
hmF2 so that they can be used in the FLIP model when measured NmF2 and hmF2 are not available. As 
a result, the FLIP model now has the capability of reproducing the IRI NmF2 and hmF2 while producing 
more physically realistic topside densities. The FLIP model is run in the normal mode without NmF 2 
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normalization when studying 'low well the current theory allows the model to reproduce the electron density 
data, while in those studies that demand accurate electron densities as ancillary information, the NmF 2 
normalization is invoked. 

As mentioned in our previous annual report, a study of the ionosphere at solar maximum was conducted 
with the original intention of ascertaining the effects of vibrationally excited N 2 (NJ) on F-region ion densities 
because our previous studies i ndicated that vibrationally excited N 2 should be important in summer at solar 
maximum but not in winter. We performed comparisons between our FLIP model and a large amount of 
Millstone Hill digisonde and Hobart ionosonde data for 1990 when the solar activity was generally very high 
but also varied by a factor of 2 We found that the measured densities were almost insensitive to solar activity 
while the model had a linear dependence on solar activity. When vibrationally excited N 2 was included in 
the calculation, the model and data were brought into agreement. 

Further investigation in the past 6 months revealed inadequacies in our EUV flux model as well as 
in other commonly used solar EUV flux models. As a result of a detailed study of available EUV flux 
data, a new solar extreme ultraviolet flux model, based on the F74113 reference spectrum and the AE-E 
relative flux measured variations, has been developed and a paper is now in press. The new EUV flux model 
(EUVAC) produces 50-575 A fluxes that are in excellent agreement with the most reliable measurements 
(Figure 2). The photoelectron fluxes produced using EUVAC are also in excellent agreement with measured 
photoelectron fluxes. The variation of the fluxes with daily F10.7 from EUVAC is smaller than the fluxes 
from the older model. This is because the long term average F10.7 flux was constant even though the daily 
F10.7 index changed a lot in 1990. With the development of the new EUV flux model there was excellent 
agreement between the modeled and measured NmF 2 at both Hobart and Millstone Hill in 1990 without 
the need to invoke N?. The observed lack of variation in NmF2 as a function of daily F10.7 is partly due 
to the relationship between F 10.7 and the solar EUV flux but is also due to changes in the thermosphere. 
X* vvas not needed to produce good agreement between model and data under these conditions and further 
study will be needed to quan lfy its role in the ionosphere. This study has been reported in a paper that is 
in press. 

An undergraduate student (Trevor Garner), has been doing research in regard to anomalous electron 
temperature enhancements over Millstone Hill. Trevor has done a detailed study of the occurrence char- 
acteristics of these events and found that it is mostly a winter phenomenon which occurs on 70% of the 
January nights and rarely in summer. However, it also occurs from September through May. The events 
typically occur prior to midright and last for 2-3 hours. This behavior seems to rule out a direct link to 
conjugate photoelectron heating which should be absent at the equinoxes and present all night in January. 
He presented a poster paper on this subject at the 1992 CEDAR meeting and also at the NASA sponsored 
Magnetospheric Models conference held at Guntersville, Alabama in October 1993. A paper submitted to the 
Journal of Geophysical Research is in press. The newly developed capability of the FLIP model to reproduce 
measured ;VmF2 and hmF 2 will be extremely useful in studying the mechanism that leads to the elevated 
temperatures at Millstone Hill. We will be able to determine the characteristics of the plasmasphenc heat 
flux that causes the enhanced temperatures. This could lead to a better understanding of plasmaspheric 
processes. 

A major study has been conducted of the ionospheric behavior at Millstone Hill and Arecibo during the 
March 1990 magnetic storms It was found that the FLIP model reproduced the electron density very well 
at Millstone Hill and at Are: ibo during the the quiet periods and during the first of two storms. During 
the second storm, which occ.ired two days after the first, the FLIP model was not capable of reproducing 
the measured electron density as well. These results have been presented at a special storm workshop at 
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Millstone Hill in March 1993, at the Baltimore AGU meeting in May 1993, and at the June 1993 CEDAR 
meeting in Boulder in 1993 ar d a journal paper is in preparation. 

Another major study was conducted in conjunction with the Lower Thermosphere Coupling Study 
(LTCS). The FLIP model electron temperature and density were compared with the radar measurements 
at Millstone Hill in December 1989 and 1991. These periods were very quiet magnetically but the solar 
acivity was high. The measured and modeled densities were very close during the day but the model 
nighttime densities were low by 50% to a factor of 2. When the N mF 2 normalization was used, the model 
temperature agreed very well with the measured temperature. These results were presented in a special 
session at the December 1993 AGU meeting in San Francisco. A paper is in preparation on this topic. 

During 1993, a graduate student (Michael Zhu) developed the framework for a new empirical thermo- 
spheric wind model based on 1 he techniques I have developed [ Richards , 1991]. This work is being continued 
by Kevin Zou. The first version of the new wind model is on schedule to be completed by the end of summer 
1994. Krishna Madakasiri is developing a new model of the peak height and density of the F2 region for use 
with the Richards [1991] technique for deriving winds from hmF2 . 

An ancillary task under his grant was to make the FLIP model available to the wider community by 
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Figure 1: This figure shows the remarkably good agreement between the measured electron 
density (dashed contours) and the modeled electron density from the FLIP model using the 
new NmF 2 technique. The measurements were made at Millstone Hill from December 8 to 
II, 1991 
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Figure 2: The measured total EUV 50-575 A irradiance as a function of F10.7. See Feng 
et ai [1989] for details of the measurements and Tobiska [1991] for SERF2 model values. 
Note that we have circled the AE-E fluxes given by Feng et al. to distinguish them from 
the AFGL rocket measurements. The solid line with asterisks is the integrated flux from the 
FLIP model. The solid line labeled HFG is a weighted least squares fit to the Hintereggcr et 
ai [1981] model output for the AE-E measurement period between mid 1977 and the end 
of 1980. The broken curve is a weighted least squares fit to the output from our new model 
(EU VAC) for the same period. 
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0 + (^P) is an important species for studies of the ionosphere and thermosphere: its emission at 
7320 A can be used as a diagnostic of the thermospheric atomic oxygen density. Unfortunately, 
there are no laboratory measurements of the O and N? reaction rates which are needed to determine 
the major sinks of 0+( 2 P). The reaction rates that are generally used were determined from 
aeronomic data by Rusch et al. (1977) but there is evidence that several important inputs that 
they used should be changed. We have recalculated the O and Nj reaction rates for 0*^( 2 P) 
using recent improvements in the solar EUV flux, cross sections, and photoelectron fluxes. For 
the standard solar EUV flux, the new N 3 reaction rate of 3.4 ±1*5 X 10 10 cm^s 1 is close to the 
value obtained by Rusch et al. (1977), but the new O reaction rate of 4.0±1.9 x 10“ 10 cm 3 *" 1 is 
about 8 times larger. These new reaction rates are derived using neutral densities, electron density, 
and solar EUV fluxes measured by Atmosphere Explorer C in 1974 during solar minimum. The 
new theoretical emission rates are in good agreement with the data for the two orbits studied by 
Rusch et al. and they are in reasonable agreement with data from five additional orbits that are 
used in this study. We have also examined the effect of uncertainties in the solar EUV flux on the 
derived reaction rates and found that 15 % uncertainties in the solar flux could cause additional 
uncertainties of up to a factor of 1.5 in the O quenching rate. 


1 . INTRODUCTION 

The 7320/7330 A doublet arising from the transition 
0 + ( 3 D - 3 P) is a very important emission for studies of 
the thermosphere. For example, Fennelly et al. [1991] and 
McDadt et al . [1991] have shown how ground-based twi- 
light measurements of this feature can be used to retrieve 
altitude profiles of atomic oxygen. However, the accuracy of 
the retrieval depends on the validity of the photochemical 
model adopted for 0 + ( 3 P). 

There are no laboratory measurements for key 0 + ( 3 P) 
reaction rates and the rates that are currently used by the 
aeronomy community were derived from Atmosphere Ex- 
plorer measurements by Rusch et ai [1977] using the best 
information that was available at that time. By combin- 
ing the solar EUV fluxes and photoionisation cross sections, 
they were able to deduce chemical reaction rates for losses 
of 0 + ( 3 P) to atomic oxygen and molecular nitrogen. Over 
the past 15 years, there have been important refinements to 
some of the input parameters that were used by Rusch et al, 
[1977]. These refinements include changes to the solar EUV 
flux, photoionization cross sections, and the importance of 
photoelectron impact ionization. R. Link (private commu- 
nication, 1990) pointed out the need to update the Rusch et 
al. model to take into account these refinements. 

Copyright 1993 by the American Geophysical Union. 

Paper number 93JA00957. 

0148-0227/93/93J A-00957 $05 .00 


Problems with the solar EUV flux were discovered by 
Richards and Torr [1984] when they compared the measured 
shape of the pbotoelectron spectrum with the expected the- 
oretical shape from the F74113 solar EUV flux spectrum of 
Hinteregger ef oi. [1977]. Another factor affecting the pho- 
toionization rate is the modification of the photoionization 
cross sections for O which are lower than those used in previ- 
ous AE work. These new cross sections have been measured 
by Samson and Pareek [1985] and Angel and Samson [1988]. 

In the work of Rusch et al . [1977], ion production by pho- 
toelectron impact was thought to be small and neglected, 
but more recent work by Richards and Torr [1988] yielded 
significant production rates that ranged from 0.14 to 0.23 of 
the photoionisation rate between the altitudes from 180 km 
to 350 km. The ratio is altitude dependent because it is a 
function of EUV attenuation, photoelectron transport, and 
whether or not the conjugate ionosphere is sunlit. 

The objective of this paper is to reevaluate the reaction 
rate coefficients for 0 + ( 3 P) with O and Na in light of the 
new information on solar fluxes, cross sections, and photo- 
electron fluxes. Therefore we have reanalyzed the two orbits 
used by Rusch et al. [1977], and five additional AE-C orbits. 

2 . THEORY 

0 + ( 3 P) is produced in the ionosphere by photoionization 
and by photoelectron impact ionization of O. It is lost in 
reactions with neutrals and electrons and by radiation. The 
major source of 0 + ( 3 P) is photoionization of O by wave- 
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